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Analytical  treatment  of  the  self-phase-modulation  of  an  ultrashort  light  pulse  is  extended  beyond  the  slowly 
varying  envelope  approximation.  The  resulting  wave  equation  is  modified  to  include  corrections  to  self-phase- 
modulation  due  to  higher-order  spatial  and  temporal  derivatives.  Analytical  solutions  are  found  in  the  limiting 
regimes  of  high  nonlinearities  and  very  short  pulses.  Our  results  reveal  features  that  can  significantly  impact 
both  pulse  shape  and  the  evolution  of  the  phase. 
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I.  INTRODUCTION 

With  laser  science  breaking  the  femtosecond  barrier,  it  has 
become  more  and  more  difficult  to  justify  the  application  of 
the  slowly  varying  amplitude  approximation  (SVEA)  to  ex¬ 
plain  new  physical  phenomena.  As  light  pulses  become  in¬ 
creasingly  shorter,  new  physical  phenomena  inaccessible  to 
the  SVEA  analysis  begin  to  play  an  important  role,  and  may 
sometimes  dominate  the  dynamics  of  ultrashort  pulses  and 
field  wave  forms.  Recent  experiments  demonstrating  the 
generation  of  attosecond  pulses  [1-3]  and  measurements 
with  these  pulses  [4,5]  raise  a  number  of  equally  challenging 
and  exciting  issues  related  to  the  propagation  and  nonlinear- 
optical  interactions  of  pulses  of  extremely  short  durations.  In 
view  of  these  breakthroughs,  many  classical  results  obtained 
in  the  realm  of  ultrafast  optics  may  have  to  be  revised  by 
extending  the  analysis  of  Maxwell  equations  beyond  the 
SVEA.  This  task  becomes  especially  urgent  when  ultrashort 
pulses  undergo  nonlinear  optical  transformations  while 
propagating  over  large  distances  as,  for  example,  in  conven¬ 
tional  [6]  or  microstructure  [7-10]  fibers,  or  long-distance 
pulse  propagation  in  the  atmosphere  [11,12]. 

Recent  numerical,  non-SVEA  studies  of  propagation  and 
nonlinear  optical  interactions  of  ultrashort  light  pulses  have 
revealed  new  interesting  scenarios  of  spatiotemporal  dynam¬ 
ics  and  spectral  evolution  of  ultrashort  pulses  [13-16],  new 
phenomena  related  to  the  ultrafast  nonlinear  response  of  a 
medium  [17,18],  and  the  influence  and  interplay  of  propaga¬ 
tion  and  phase-matching  effects  [19],  to  name  a  few.  In  the 
absence  of  analytical  solutions  it  therefore  seems  appropriate 
and  timely  to  attempt  to  put  some  of  these  results  on  a  firmer 
theoretical  ground  with  a  systematic  analytical  treatment  of 
non-SVEA  effects.  Accomplishing  this  task  involves  the  re¬ 
vision  of  the  equations  of  motion  for  basic  nonlinear  optical 
processes,  such  as  self-  and  cross-phase  modulation,  har¬ 
monic  generation,  four-wave  mixing,  and  others,  with  the 
inclusion  arising  from  the  inclusion  of  second-order  longitu¬ 
dinal,  spatial,  and  temporal  derivatives. 

In  this  paper,  we  provide  the  non-SVEA  analysis  specifi¬ 
cally  for  self-phase-modulation  (SPM)  as  a  part  of  this  plan, 
and  find  analytical  solutions  for  the  limiting  regimes  of  high 
nonlinearities  and  very  short  pulses.  Close  scrutiny  of  these 


solutions  then  reveals  important  features  in  the  dynamics  of 
the  pulse  envelope  and  phase  related  to  non-SVEA  effects. 


II.  REVISED  EQUATION  FOR 
SELF-PHASE-MODULATION 


We  begin  our  analysis  by  writing  the  scalar  wave  equation 
for  the  electric  field  E  in  an  isotropic  medium  with  no  dis¬ 
persion, 


n1  cpE  4'rrd1PIIi 

~2  n.2  =  ~2  n.2 

at  c  at 


(1) 


where  n  is  the  refractive  index  and  c  is  the  speed  of  light  in 
the  vacuum. 

Without  loss  of  generality,  a  short  pulse  propagating  along 
the  z  axis,  of  carrier  frequency  co  such  that  k=am/c,  can  be 
expressed  as  follows; 

E(x,y,z,t)  =  s(x,y,z,t)eiikz-M,)  +  c.c.  (2) 

We  write  the  nonlinear  polarization  of  the  medium  in  the 
form  of  a  power-series  expansion,  keeping  the  terms  up  to 
the  fifth  order  in  the  field  E: 

Pnl  =  X0)\E\2E  +  X(5)\E\4E,  (3) 

where  x^  and  a'5*  are  the  third-  and  fifth-order  nonlinear- 
optical  susceptibilities,  respectively.  In  writing  Eq.  (2),  we 
assume  that  the  carrier  frequency  u>  can  be  defined  for  our 
light  pulse,  which  is  true  even  for  few-cycle  pulses.  In  fact, 
Eq.  (2)  is  quite  general  in  that  it  is  equivalent  to  a  coordinate 
transformation.  Expression  (3)  corresponds  to  the  usual  re¬ 
gime  of  perturbative  nonlinear  optics.  However,  our  analysis 
will  go  beyond  the  standard  SVEA  treatment  of  self-phase- 
modulation.  We  will  find  non-SVEA  corrections  to  the  SPM- 
induced  nonlinear  phase  shift  and  show  that  these  corrections 
can  be  described  in  terms  of  the  effective  fifth-order  nonlin¬ 
ear  susceptibility.  Since  the  term  with  the  “true”  fifth-order 
susceptibility  is  kept  in  the  expression  for  the  nonlinear  po¬ 
larization  [Eq.  (3)],  we  will  be  able  to  compare  ^5^-like  non- 
SVEA  corrections  with  the  contribution  of  the  true  fifth- 
order  nonlinearity. 

Substituting  Eqs.  (2)  and  (3)  into  Eq.  (1)  for  a  plane  wave, 
we  arrive  at 
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2  cPe  de  n2  cPe  n2  d  e 

V  T  e  +  — ~  +  2 ik —  -  — =■ — ~  +  2 iu>—  — 
dz2  dz  c2dt2  c2  at 


47 T 
~JX 


,(3)I 


d2  d  ,  I  ,  4tt 

|s|e+“* 


.(5) 


X)  — r  -  2/w —  -  or  f  ere. 


df 


dt 


(4) 


To  represent  Eq.  (4)  in  dimensionless  form,  we  introduce 
the  nonlinear  length  l/Ln/=(47r<A>2|e0|2/2kc2)y3\  the  nor¬ 
malized  running  time  r=(t-nz/c)/ t0  (where  r0  is  the  pulse 
duration),  and  the  dimensionless  propagation  coordinate  £ 
=  z!Lnl.  We  normalize  the  field  with  respect  to  its  maximum 
amplitude  (e— >e/e0).  With  these  transformations  and  scal¬ 
ings  and  with  an  assumption  that  transverse  effects  are  neg¬ 
ligible,  i.e.,  V2±e  =  0,  Eq.  (4)  can  be  rewritten  as 

de  1  [  X  ^e  X  gs  } 

'di;  +  4tt[l„/(?£2  crQd  rd  £;  J 

[  \ 27T/  \ctq)  dr-  ttct0  St  J 

|gn|2y(5) 

..i  I?  I  01  A  |  id 

Xe  e- - — e  e.  (5) 


Now  we  introduce  the  parameters  fx:=\l\'nlXILnl  and 
/Ar=  1/4t7-/A./c70,  representing  the  length  and  time  scales,  re¬ 
spectively.  Under  normal  circumstances,  the  nonlinear  length 
is  much  larger  than  the  wavelength  of  light.  By  the  same 
token,  even  for  a  pulse  only  a  few  optical  cycles  in  duration, 
//-  can  be  much  smaller  than  unity.  Therefore,  keeping  only 
terms  up  to  the  first  order  in  /a,  and  /xT  in  Eq.  (5),  we  derive 

de  c^e  J  c^e  f?(|e|2e)|  , 

i —  =  -  ix- — 2  +  ixT\ - -  4 i - r  -  she 

di;  dg  [dijdT  dr  J 

2y(5) 

7^— |e|4e.  (6) 


We  now  focus  our  attention  on  the  higher-order  spatial 
and  spatiotemporal  derivatives  that  appear  on  the  right-hand 
side  of  Eq.  (6).  Differentiating  Eq.  (6)  with  respect  to  z,  and 
neglecting  the  terms  of  the  order  of  /x2,  /a2,  /xjx-,  and 
we  find 


cPb 


-  d2 s  .  .  <?|e|  7de 

T^.  =  ltlz~Ts  +  l2\s\s^7  +  l\s\  —t  +  0{fiT) 
d  f  d  %  d  E,  d  f 

=  -  |e|4e  +  0(/xT)  +  0(fxz) . 


(7) 


To  estimate  d\s\l di;,  we  represent  the  field  as  s=|e|e'^  and 
substitute  this  expression  into  Eq.  (6),  which  yields  d\s\l di; 
=  0(/xJ  +  0(/x,).  To  find  the  derivative  d2B,l di; dr,  we  once 
again  differentiate  Eq.  (6),  this  time  with  respect  to  r, 


c?s 
d  £;  d  T 


+  0(fxT)  +  0(/ttz)]. 


(8) 


Then,  substituting  Eqs.  (7)  and  (8)  into  Eq.  (6),  and  keeping 
the  terms  of  the  zeroth  and  first  orders  in  /xT  and  /a,,  we  find 


de  4 

i —  =  a.  e  e  - 

di  1 


yv5),  |4 

- At—  e  s 

X(> 


1  +  3iuT —  r  sre. 
OtV 


(9) 


Expression  (9)  gives  a  SPM  equation  revised  to  include 
non-SVEA  corrections  related  to  higher-order  derivatives  in 
Eq.  (1),  as  well  as  the  term  related  to  the  fifth-order  nonlin¬ 
earity.  With  /x-=/xT=0,  Eq.  (9)  reduces  to  the  standard, 
SVEA  equation  for  self-phase-modulation.  The  revised  SPM 
equation  (9)  is  instructive  from  the  physical  point  of  view,  as 
it  provides  useful  insights  into  the  physical  origin  of  the 
first-order  non-SVEA  corrections  and  allows  these  correc¬ 
tions  to  be  compared  with  the  contribution  of  the  fifth-order 
nonlinearity.  The  term  involving  jxz  gives  a  ^5)-like,  i.e., 
quintic  in  the  field,  correction  to  the  SPM  equation.  This 
finding  is  fundamental  from  the  viewpoint  of  the  basic  prin¬ 
ciples  of  nonlinear  optics,  as  it  reveals  the  existence  of  an 
additional  physical  channel  whereby  effects  quintic  in  the 
field  may  contribute  to  the  SPM  of  ultrashort  laser  pulses.  In 
the  presence  of  diffraction,  for  example,  a  ^<5)- like  term 
would  tend  to  defocus  the  beam,  leading  to  increased  non¬ 
linear  thresholds  and  modified  spatiotemporal  dynamics  for 
long  distance  propagation  [15].  It  can  be  easily  seen  from  the 
definition  of  fx;  that  the  role  of  the  ^,5l-lil<c  term  in  Eq.  (9) 
increases  with  the  decrease  in  the  nonlinear  length  Lnh  i.e., 
with  the  growth  in  the  nonlinearity.  Thus,  y5,-|i|<e  effects 
may  influence  nonlinear-optical  interactions  governed  by  the 
third-order  nonlinearity  even  in  the  regimes  of  perturbative 
nonlinear  optics.  The  term  involving  /at  gives  a  small  correc¬ 
tion  to  the  term  proportional  to  |e|2e,  which  dominates  the 
SPM  in  the  SVEA  regime.  This  correction,  as  it  follows  from 
the  definition  of  the  parameter  /x„  becomes  significant  in  the 
case  of  very  short,  few-  and  single-cycle  pulses.  As  will  be 
shown  in  the  following  sections,  this  correction  may  also 
give  rise  to  noticeable  changes  in  the  pulse  envelope  and 
phase  distribution  evolution  of  self-phase-modulated  pulses. 

Expression  (9)  involves  two  terms  quintic  in  the  field  E. 
The  physical  origin  of  these  two  terms  is  different.  While  the 
first  term  on  the  right-hand  side  originates  from  the  non- 
SVEA  correction  related  to  the  second  derivative  in  z,  the 
second  term  describes  the  contribution  of  the  true  fifth-order 
nonlinearity  of  the  medium.  It  would  be  instructive  to  gain  a 
general  understanding  of  the  relative  significance  of  these 
two  terms  in  Eq.  (9)  by  using  the  relation  between  off- 
resonance  nonlinear  susceptibilities  of  different  orders  typi¬ 
cal  of  the  perturbative  regime  of  nonlinear  optics  [20]: 
X(n)~XW/\Eat\n-\  where  Eat  is  the  characteristic  atomic  field 
(typically  on  the  order  of  108— 109  V/cm).  The  first 
two  terms  on  the  right-hand  side  of  Eq.  (9)  can  be 
then  estimated  as  leoiy5^Y(3)~leo|2/|£'ar|2  and 
/Az(7r|80|2/n2)y3)~(7r/n2)yi)(|80|2/|Ej2)~|80|2/|£j2.  The 
non-SVEA  correction  in  a  nonresonant  situation  is  thus  of 
the  same  order  of  smallness  as  the  term  related  to  the  true 
fifth-order  nonlinearity.  These  two  channels  of  nonlinear- 
optical  interactions  can,  therefore,  interfere  either  construc¬ 
tively  or  destructively,  giving  rise  to  interesting  effects,  en¬ 
hancing  or  suppressing  nonlinear  phenomena,  depending  on 
the  relative  sign  of  the  first  two  terms  in  Eq.  (9).  In  what 
follows,  we  will  present  an  iterative  procedure,  allowing  an 
analytical  integration  of  non-SVEA  SPM  equations  for  the 
amplitude  and  the  phase  of  ultrashort  pulses  including  only 
the  cubic  optical  nonlinearity. 


053803-2 


ANALYTICAL  TREATMENT  OF  SELF-PHASE-... 


PHYSICAL  REVIEW  A  69,  053803  (2004) 


III.  ANALYTICAL  SOLUTIONS  TO  THE  NON-SVEA 
SPM  EQUATION 

In  this  section,  we  will  develop,  following  the  plan  out¬ 
lined  above,  an  iterative  procedure  solving  the  generic  non- 
SVEA  SPM  equation  (9)  for  a  medium  with  cubic  and  quin- 
tic  nonlinearities.  Representing  the  field  of  a  light  pulse  as 
e(£,  r)=A(£,  r)exp[/0(£,  r)],  we  arrive  at  the  following  set  of 
equations  for  the  amplitude  A  and  the  phase  <f>  of  the  field: 


dA 

H 


=  -9/jltA2 


dA 
d  t’ 


phase.  Substituting  Eq.  (16)  into  the  second  equation  in  the 
set  (10),  we  derive  the  following  equation  for  the  first-order 
non-SVEA  correction  to  the  phase: 


^(o  =  _  2A?d<ho) 

dg  dr 


(18) 


Substituting  the  first-order  non-SVEA  solution  (15)  for 
the  amplitude  into  Eq.  (18),  keeping  only  the  terms  of  the 
zeroth  and  first  orders  in  /jlt  and  jl.  in  the  resulting  equation, 
and  assuming  that  the  input  pulse  has  no  initial  chirp,  we 
arrive  at  the  following  solution  for  the  phase: 


d  <fi 


=  A2-jH-A4-3ii7A2—, 

d  t 


where  ^z= ^z-\eQ\2 X(S) I X0) ■ 

We  now  write  the  field  amplitude  as  a  series  expansion  in 
the  small  parameter 


A -A(o)  +  nA(i)  +  +  ' ' '  >  (11) 


where  A(0)  is  the  SVEA  solution  to  the  SPM  equation  for  the 
field  amplitude. 


dA(0) 


(12) 


and  A(!)  and  Al2)  are  the  non-SVEA  corrections  of  the  first 
and  second  orders  in  respectively.  The  equation  for  the 
first-order  non-SVEA  correction  to  the  field  amplitude  is 
then  written  as 


^0)(i,T)  =  [A2-jd7A40-\H/x^ 


l3dAQ(r) 


d  t 


0u)(£  t)  =  -  6Ao 


dA0(r) 
d  T 


e. 


Combining  the  solutions  for  the  phase  and  the  amplitude 
given  by  Eqs.  (15)  and  (19),  we  derive  the  following  expres¬ 
sion  for  the  field: 


e(£,r)  =  A0 


1  —  9/x^40 


^4q 

d  T 


exp[«'(0svEA  -  A<k  -  A (f>T)], 


(20) 


where 


^svea-Ao^.  (21) 

is  the  SVEA  solution  for  the  SPM-induced  phase  shift. 


aAW 


=  -9  A(0)2 


^A(o) 
d  T 


(13) 


The  solution  to  Eq.  (12)  is  well  known  [20]:  A(0)(£,  r) 
=A0(t),  where  A0(r)  is  the  initial  pulse  envelope  normalized 
to  the  maximum  field  amplitude  [A0(r)  =  exp(-T2)  for  a 
Gaussian  pulse].  Integration  of  Eq.  (13)  yields 


A(1)(£t)  =  -9A5-^.  (14) 

d  T 


The  solution  to  the  non-SVEA  SPM  equation  (9)  in  the  first- 
order  approximation  in  juT  is  thus  written  as 

1  d  Af| 

A(£7-)=A„-9/V^— A0  =  A0(0,  r) .  (15) 

d  T 


We  can  now  proceed  with  the  solution  of  the  equation  for 
the  phase  in  the  set  of  Eqs.  (10).  We  represent  the  phase  as 

(f>=  cf>(0)  + t*T<f>(i)  + fST(f>(2)+ ■■■  (16) 

and  take  the  solution  to  the  equation 


n 


=  A2-/HA4 


(17) 


as  the  zeroth-order  iteration.  In  view  of  Eq.  (17),  the  small 
parameter  Jiz  also  appears  in  the  series  expansion  (16)  for  the 


A^=24  (22) 

d  t 

is  the  non-SVEA  correction  to  the  nonlinear  phase  shift  re¬ 
lated  to  the  second  derivative  in  time,  and 

A  </>.=  jiAte  (23) 

is  the  non-SVEA  correction  originating  from  the  second  de¬ 

rivative  in  z. 

For  a  pulse  with  an  initial  phase  modulation  <p(0 ,  r) 
=  </>o(t)  =  ar  (afiT<t  1 ),  the  correction  to  the  nonlinear 
phase  shift  is  given  by 

A<f>T=24piAl—^e-6tJLA2Qa.  (24) 

d  T 

Expression  (23)  suggests  that  the  non-SVEA  correction 
A  cj).  to  the  SPM-induced  phase  shift  is  quartic  in  the  field. 
This  functional  dependence  is  typical  of  nonlinear- 

optical  effects.  In  our  case,  however,  the  smallness  of  the 
phase  shift  A  <pz  is  related  not  only  to  the  fifth-order  suscep¬ 
tibility  x^5\  but  is  governed  by  the  parameter  jlz,  which  is  in 
turn  controlled  by  both  the  third-  and  fifth-order  susceptibili¬ 
ties  ^3)  and  y'51.  In  a  practically  important  situation  of  a 
medium  with  a  resonance-enhanced  third-order  susceptibility 
X®,  the  phase  shift  A<b:,  given  by  Eq.  (23)  can,  therefore, 
noticeably  exceed  the  true  y^-induccd  phase  shifts,  giving 
rise,  as  will  be  shown  below,  to  detectable  (and  sometimes 
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FIG.  1.  The  nonlinear  phase  shift  </>SVEA  calculated  within  the 
framework  of  the  slowly  varying  envelope  approximation  (solid 
line)  and  the  non-SVEA  correction  A </>.  to  the  nonlinear  phase  shift 
related  to  the  second-order  derivative  in  z  (dashed  line)  for  £  =  1 
with  /Z-=10-3.  The  variable  t  is  normalized  to  the  pulse  duration. 

significant)  ^5,-like  corrections  to  the  phase  in  the  purely 
perturbative  regime  of  nonlinear-optical  interactions. 

Non-SVEA  effects  related  to  the  second-order  time  de¬ 
rivative,  as  can  be  seen  from  Eqs.  (22)  and  (24),  give  rise  to 
distortions  in  the  envelope  of  a  light  pulse,  which  becomes  z 
dependent  even  in  the  absence  of  group-velocity  dispersion. 
Expressions  (22)-(24)  are  also  instructive  in  demonstrating 
that  both  pulse-envelope  and  phase-profile  distortions  related 
to  the  non-SVEA  effects  accumulate  with  z  as  the  pulse 
propagates  through  the  nonlinear  medium.  This  dependence 
on  the  longitudinal  coordinate  can  lead  to  noticeable  devia¬ 
tions  from  the  SVEA  results  for  pulse  envelopes  and  nonlin¬ 
ear  phase  shifts  for  sufficiently  large  propagation  distances, 
as,  for  example,  in  the  case  of  standard  [6]  or  microstructure 
optical  fibers  [7-10]  or  nonlinear-optical  interactions  in  the 
atmosphere  [11,12]. 

IV.  RESULTS  AND  DISCUSSION 

To  assess  the  role  of  non-SVEA  effects  in  self-phase- 
modulation,  we  use  the  solutions  to  the  revised  SPM  equa¬ 
tion  derived  above  to  calculate  the  non-SVEA  phase  shifts 
and  distortions  in  pulse  envelopes  for  different  values  of  pa¬ 
rameters  Jlz  and  ilt.  Figure  1  displays  the  nonlinear  phase 
shift  </>svea  calculated  within  the  framework  of  the  standard 
SVEA  approach  (solid  line)  and  the  non-SVEA  correction 
A  <!>,.  to  the  nonlinear  phase  shift  related  to  the  second-order 
derivative  in  z  (dashed  line)  for  t=  1 .  The  non-SVEA  phase 
shift  was  calculated  using  Eq.  (23)  with  jxz  ~  1  0-3  (which 
corresponds  to  Ln/=  100  fim  in  a  medium  with  ,y(:,  =  0).  Such 
values  of  Jiz  can  easily  be  achieved,  for  example,  by  propa¬ 
gating  light  pulses  with  a  wavelength  of  1  jam  and  the  inten¬ 
sity  on  the  order  of  1012  W/cm2  through  fused  silica  (in  the 
absence  of  true  fifth-order  nonlinearity).  The  non-SVEA  cor¬ 
rection  to  the  phase  shift,  A 4>z,  as  can  be  seen  from  Fig.  1, 
can  reach  a  maximum  of  about  0.1  %  under  these  conditions. 
This  part  of  the  phase  shift,  however,  increases  proportion¬ 
ally  to  jaz,  reaching  0.1  %  for  /I. =0.01 .  These  estimates 


FIG.  2.  The  envelope  of  a  self-phase-modulated  pulse  calculated 
within  the  framework  of  the  slowly  varying  envelope  approxima¬ 
tion  (solid  line)  and  with  the  use  of  Eq.  (20)  (dashed  line),  including 
effects  related  to  the  second-order  derivative  in  r,  for  £=1  with 
jit=10~2,  jlz> jjlt.  The  variable  t  is  normalized  to  the  pulse 
duration. 

show  that  non-SVEA  phase  shifts  may  play  an  important 
role,  for  example,  in  nonlinear-optical  experiments  in  micro¬ 
structure  fibers  [10],  where  light  intensities  exceeding 
10°  W/cm2  are  often  achieved  by  coupling  femtosecond 
pulses  in  a  few-micron  fiber  core. 

In  the  regime  of  very  short  pulses,  with  jaz  8>  jxr  non- 
SVEA  pulse-envelope  and  phase-profile  distortions  accumu¬ 
late  as  the  pulse  propagates  through  a  medium.  Figure  2 
shows  distortions  in  the  envelope  of  a  light  pulse  with  an 
initial  pulse  duration  of  about  30  fs  caused  by  effects  related 
to  the  second-order  time  derivative  in  the  wave  equation.  The 
changes  in  the  pulse  envelope  are  already  quite  noticeable 
for  £=1.  The  influence  of  non-SVEA  phase  shifts  is  illus¬ 
trated  in  Fig.  3,  which  compares  the  nonlinear  phase  shift  of 
a  pulse  with  a  duration  of  about  30  fs  calculated  within  the 
framework  of  the  SVEA  approach  with  the  nonlinear  phase 
shift  calculated  with  the  use  of  Eqs.  (20)  and  (22),  including 


T 


FIG.  3.  The  nonlinear  phase  shift  </>SVEA  calculated  within  the 
framework  of  the  slowly  varying  envelope  approximation  (solid 
line)  and  the  correction  to  the  phase  shift  A  cf>T  calculated  with  the 
use  of  Eq.  (22)  (dashed  line),  including  effects  related  to  the 
second-order  derivative  in  r,  for  £=1  with  jjlt=  10-2.  The  variable  t 
is  normalized  to  the  pulse  duration. 
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FIG.  4.  The  nonlinear  phase  shift  of  a  short  pulse  with  an  initial 
linear  chirp  calculated  within  the  framework  of  the  slowly  varying 
envelope  approximation  (solid  line)  and  the  correction  to  the  phase 
shift  A(f>T  calculated  with  the  use  of  Eq.  (24)  (dashed  line),  includ¬ 
ing  effects  related  to  the  second-order  derivative  in  t,  for  £=  1  with 
o:=-0.5  and  fiT=  10-2. 

non-SVEA  corrections  related  to  the  second-order  time  de¬ 
rivative.  Similar  to  pulse-envelope  distortions,  the  deviation 
of  the  non-SVEA  nonlinear  phase  shift  from  the  SVEA  result 
becomes  quite  noticeable  already  for  £=1,  increasing  in  the 
process  of  pulse  propagation  in  accordance  with  Eqs.  (20) 
and  (22). 

Putting  an  initial  chirp  on  the  pulse  adds  phase-control 
aspects  to  the  analysis  of  nonlinear  interactions  of  ultrashort 
pulses.  Figures  4  and  5  show  that  such  an  analysis  should 
include  non-SVEA  effects,  which  may  have  a  noticeable  in¬ 
fluence  on  phase-controlled  nonlinear-processes.  Figure  4 
compares  the  SVEA  nonlinear  phase  shift  of  a  short  pulse 
having  an  initial  linear  chirp  with  a  chirp  parameter  a= 
-0.5  and  //_=  1  0~2  with  the  results  of  calculations  performed 
with  the  use  of  Eq.  (24),  including  effects  related  to  the 
second-order  time  derivative.  The  evolution  of  the  nonlinear 
phase  shifts  for  transform-limited  and  initially  chirped  pulses 
illustrated  in  Figs.  5(a)-5(c)  shows  how  the  non-SVEA  part 
of  the  nonlinear  shift  is  accumulated  and  the  distribution  of 
the  nonlinear  shift  within  the  pulse  becomes  more  and  more 
asymmetric  as  the  pulse  propagates  through  the  nonlinear 
medium.  Non-SVEA  phase  corrections  should  be,  therefore, 
taken  into  consideration  in  coherent-control  experiments 
[21-26],  where  chirped  ultrashort  laser  pulses  are  used  to 
steer  molecular  excitations,  wave-packet  motions,  chemical 
reactions,  or  nonlinear-optical  interactions.  Non-SVEA  ef¬ 
fects,  as  can  be  seen  from  Fig.  5,  also  noticeably  influence 
the  SPM-phase-shift  precompensation  function  of  the  initial 
chirp,  as  the  result  of  such  a  precompensation  may  deviate 
from  expectations  based  on  the  SVEA  analysis.  Positive  ini¬ 
tial  chirp,  as  can  be  seen  from  the  comparison  of  the  evolu¬ 
tions  of  the  corrections  A<pT  for  transform-limited  and  ini¬ 
tially  chirped  pulses  shown  by  the  dashed  and  dotted  lines  in 
Figs.  5(a)-5(c),  respectively,  allows  phase-profile  distortions 
related  to  non-SVEA  effects  to  be  reduced  or  precompen¬ 
sated.  This  precompensation  is  effective,  however,  only 
within  propagation  lengths  on  the  order  of  the  nonlinear 
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FIG.  5.  The  nonlinear  phase  shift  for  an  initially  transform- 
limited  pulse  calculated  within  the  framework  of  the  slowly  varying 
envelope  approximation  (solid  lines)  and  the  correction  to  the  phase 
shift  A<f>T  calculated  with  the  use  of  Eq.  (24)  for  an  initially 
transform-limited  pulse  with  a= 0  (dashed  lines)  and  a  pulse  having 
a  linear  initial  chirp  with  a=0.5  (dotted  lines)  for  £=  1  (a),  £=2  (b), 
and  f=5  (c);  /rT=  10-2. 

length  [Fig.  5(a)].  Non-SVEA  phase-profile  distortions  accu¬ 
mulate  as  the  pulse  propagates  through  the  medium  [Figs. 
5(b)  and  5(c)],  the  first  term  in  Eq.  (24)  starts  to  dominate, 
and  the  expression  for  the  non-SVEA  phase-shift  correction 
A(f>Tis  reduced  to  Eq.  (22).  This  analysis  shows  that  a  prop¬ 
erly  chosen  initial  chirp  of  short  pulses  may  serve  as  an 
important  knob  in  phase-controlled  molecular  motions  and 
nonlinear  interactions. 
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V.  CONCLUSION 

Thus,  analysis  performed  in  this  paper  shows  that  the  sce¬ 
narios  of  nonlinear-optical  interactions  of  ultrashort  pulses 
may  sometimes  noticeably  differ  from  predictions  of  the 
slowly  varying  envelope  approximation.  We  have  derived  a 
revised  equation  for  self-phase  modulation,  including  correc¬ 
tions  related  to  higher-order  derivatives  in  the  wave  equa¬ 
tion.  We  have  also  developed  an  iterative  procedure  allowing 
the  analytical  solution  of  this  equation  for  the  limiting  re¬ 
gimes  of  high  nonlinearities  and  very  short  pulses.  Analysis 
of  these  solutions  reveals  interesting  features  in  the  pulse- 
envelope  and  phase-profile  evolution  related  to  non-SVEA 
effects  in  self-phase-modulation.  It  was  shown,  in  particular, 
that  y5)-like  effects,  i.e.,  effects  of  the  fifth  order  in  the  laser 
field,  may  influence  nonlinear-optical  interactions  governed 
by  the  third-order  nonlinearity  even  in  the  regimes  of  pertur¬ 
bative  nonlinear  optics.  Many  of  the  non-SVEA  effects  ac¬ 
cumulate  with  the  propagation  coordinate,  leading  to  notice¬ 
able  changes  in  the  pulse  envelope  and  the  phase  distribution 


in  the  laser  pulse  for  large  propagation  lengths,  e.g.,  in  opti¬ 
cal  fibers  or  in  nonlinear-optical  interactions  in  the  atmo¬ 
sphere.  Non-SVEA  phase  corrections  and  pulse-envelope 
distortions  should  be  taken  into  consideration  in  coherent- 
control  experiments,  where  chirped  ultrashort  laser  pulses  are 
used  to  steer  molecular  excitations,  wave-packet  motions, 
chemical  reactions,  or  nonlinear-optical  interactions. 
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